Chain-reaction cascades in surfactant monolayer buckling 
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Certain surfactant monolayers at the water-air interface have been found to undergo, at a critical 
surface pressure, a dynamic instability involving multiple long folds of micron width. We exploit the 
sharp monolayer translations accompanying folding events to acquire, using a combination of fluo- 
rescence microscopy and digital image analysis, detailed statistics concerning the folding dynamics. 
The motions have a broad distribution of magnitudes and narrow, non-Gaussian distributions of 
angles and durations. The statistics are consistent with the occurrence of cooperative cascades of 
folds, implying an autocatalytic process uncommon in the context of mechanical instability. 

PACS numbers: 68.18.Jk, 64.60.Qb, 82.60.Nh, 87.68.-|-z 
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Surfactant monolayers are found in many systems con- 
taining water-air or water-oil interfaces where surface 
tension, wetting, or liquid-film stability are to be con- 
trolled m . In recent years a remarkable variety of three- 
dimensional structures have been discovered upon lateral 
compression of surfactant monolayers, including straight 
folds [^H,^, convoluted foldsHl^, and attached vesic- 
ular objects of various shapes 'Sj. These instabilities have 
distinctive length scales ranging between 0.1 and 10 /im. 
Thus the predominant stress relaxation is neither at the 
molecular level, which would lead to breakage or disso- 
lution 6], nor at the macroscopic one, which would lead 
to long- wavelength buckling [2|. Moreover, the transi- 
tions occur under a net tensile stress, i.e., at surface pres- 
sures smaller than the bare surface tension of water (72 
mN/m). Though the actual relaxation mechanism is im- 
known, a plausible driving force may be bilayer cohesion, 
i.e., the preference of the hydrophobic surfactant tails to 
join rather than remain in contact with air 0. Domain 
boundaries in biphasic monolayers or grain boundaries 
in monophasic ones entail nanoscale topographies, which 
should cause localized buckling at low enough, positive 
tension 0|, thus lowering the nuclcation barrier for bi- 
layer cohesion. 

Out of the newly observed structures, the straight 
folds stand out as essentially different, corresponding to 
a more solid-like, nonequilibrium response of the mono- 
layer H, H, ^ Hi ' They are observed at lower tempera- 
ture and higher compression rate, occur at higher critical 
pressures 70 mN/m), and are anisotropic, i.e., aligned 
on average perpendicular to the compression direction. 
A fold comprises a piece of monolayer of micron width 
and macroscopic length probably bound into a bilayer 
strip (Fig. The main obstacle in studying fold for- 
mation is that it is a nucleated event initiated at an un- 
predictable spot and lasting a fraction of a second. Very 
rarely is such an event captured inside the field of view 
as in Fig. ^ Nevertheless, whenever a fold forms else- 
where, the viewed piece of monolayer translates sharply 



and uniformly. Watching the monolayer jump as a result 
of such events, one is struck by the uniformity and uni- 
directionality of the motion, in contrast with the mono- 
layer heterogeneity. The statistical analysis presented be- 
low corroborates this impression, revealing such features 
as anomalously narrow distributions of translation angle 
and duration. 
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FIG. 1: (a)-(c) Fluorescence micrographs separated by 1/30 
s intervals, showing the nearly- simultaneous formation of two 
folds. The images are blurred by monolayer motion. The 
scale bar length is 50 ^m. (d) Typical output of the tracking 
program, showing the monolayer translation within the field 
of view in sequential video frames. The spikes correspond to 
folding events occurring out of view. The dotted line shows 
the threshold used for event identification. 



Insoluble monolayers containing a 7:3 mixture 
of dipalmitoylphosphocholine (DPPC) and palmitoy- 
lolcoylphosphoglycerol (POPG) were spread on the air- 
water interface from a chloroform solution at 25° C. This 
was done on a rectangular Langmuir trough of maximum 
area 145 cm^, fitted with two mobile Teflon barriers of 
length 6.35 cm each, allowing for symmetric lateral com- 
pression. The phospholipids formed a low-density mono- 
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layer (area per molecule a > 100 A^), exhibiting coexis- 
tence of two-dimensional gas and liquid phases. Lateral 
compression led, at a ~ 67 A^, to a homogeneous hq- 
uid phase followed, for a ^ 65 A^, by nucleation and 
growth of flower-like condensed domains (Fig. ^) until, 
at a ~ 30 A^, folding began. The relative barrier ve- 
locity was 0.1 mm/s, corresponding to compression rate 
of 6.35 mm^/s and strain rate at the onset of folding of 
0.00154 s~^. Throughout the compression the morphol- 
ogy was observed using cpifluorescence video microscopy 
and the surface tension monitored by a Wilhelmy surface 
balance. Further details of the apparatus can be found 
in Ref. 0. The analog video was digitized into a series 
of 8-bit grayscale bitmaps of 640 x 480 pixels at a rate of 
29.97 frames/s. The series of images were then analyzed 
using a custom-made tracking program whose algorithm 
will be detailed elsewhere A typical output is shown 
in Fig. njd). Using a velocity threshold of 2 pixels per 
frame we identified 1817 events, recording for each its 
starting time, duration and total I = (Ix, ly) translation. 

Figure |21[a) shows the distribution of translations, ex- 
hibiting a broad tail — translations ten times the most 
probable one were observed. The distribution of off-times 
t (waiting times between events) is presented in Fig.[2Ib). 
The mean off-time is (t) = 0.31 ± 0.01 s, i.e., there are 
about three events per second. The histogram fits well 
an exponential distribution, pt{t) = (i)^^e^*/^*^, consis- 
tent with a Poissonian, uncorrelated sequence of events. 
This conclusion is strengthened by a lack of correlation 
between I and the off-times before or after the event. 

Figure Ola) shows the distribution of translation an- 
gles \6\, where tan 61 = ly/lx, and 6 = corresponds to 
motion parallel to the compression direction. The distri- 
bution is sharply peaked at = 0, with standard devia- 
tion ae = {9^y/^ = 16.0 ± 0.3°. Folding is thus highly 
anisotropic, implying an elastic response of the mono- 
layer within the folding time scale. This is in line with 
the viscoelasticity revealed bysurface-rheology measure- 
ments in similar systems p|, llll , which yields relaxation 
times of order tens of seconds. One would expect a Gaus- 
sian distribution of fold angles — because of either the 
combined effect of many scattering factors or a Boltz- 
mann factor for the nucleation of a slanted fold with 
respect to the direction of maximum stress, whose en- 
ergy increases as sin^ 9^6^. The peak of the measured 
distribution, however, is much sharper and clearly can- 
not be fitted by a Gaussian distribution (dotted curve). 
In Fig. inib) we show the mean angle \9i\ of events as a 
function of their translation /. The two quantities are 
anticorrelated — larger events have smaller angles. 

Figure Ela) presents the distribution of on-timcs T 
(event durations). The distribution is narrow and asym- 
metric, yielding a mean on-time (T) = 0.124 ± 0.001 s 
and standard deviation ctt — 0.051 ± 0.001 s. The corre- 
sponding Gaussian distribution is depicted by the dotted 
line, highlighting the anomalous shape of the measured 
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FIG. 2: (a) Distribution of translations. The solid line is 
obtained assuming a cascade mechanism [Eq. Q]. (b) Dis- 
tribution of off-times. The solid line shows an exponential 
distribution using the measured mean off-time (t) = 0.31 s. 



distribution. The more moderate decrease to the right of 
the peak fits an exponential decay rather than a Gaus- 
sian one (inset). The narrow T distribution is surprising 
in view of the broad / distribution; one expects larger 
events to last longer. This correspondence is verified in 
Fig. Hb) , where the average on-time T/ is plotted as a 
function of the translation I. Yet, the increase of Ti with 
I is only logarithmic and, therefore, even very large trans- 
lations do not have correspondingly long duration. 

The broad distribution of translations can arise from 
two alternative scenarios: each observed event could cor- 
respond to either a single fold, whereby the folds have a 
broad distribution of sizes, or a cascade of roughly iden- 
tical folds, the cascades having a broad distribution of 
magnitudes. The following analysis, as well as the hand- 
ful of folds captured in the field of view, strongly support 
the latter scenario. 

Let us assume that an observed event is caused by a 
cascade of n folds, each contributing roughly the same 
translation li and having an angle 9i, i = I . . .n. We 
assume for simplicity that the angles are drawn from in- 
dependent Gaussian distributions having a standard de- 
viation aei, VeM) = (27rCT^i)-i/2g-e,V(2<T,\)^ rpj^g 

sured translation is the sum of contributions from all 
folds, / ~ nil, but the observed translation angle is the 
average of the individual fold angles, 9 ~ "'~^X]r=i^«- 
The distribution of 9 is therefore narrower the larger the 
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FIG. 3: (a) Distribution of translation angles. The dotted line 
shows a Gaussian distribution with the measured (S^). The 
solid line is obtained assuming a cascade mechanism [Eq. ||3J] 
with aei = 16°. (b) Average angle of events whose transla- 
tions fall within the same narrow range. The solid line is a fit 
to \ei\ = 7ri/2 [Eq. Q] with 7 = 18 ± 3 Ofim^^'^. The in- 
set presents the same data on a log-log scale, the dashed line 
being an error- weighted linear fit with a slope of —0.44 ±0.05. 
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FIG. 4: (a) Distribution of on-times. The dotted line shows a 
Gaussian distribution using the measured mean and standard 
deviation. The inset presents the data on a linear-ln scale, 
the solid line being an error-weighted linear fit with a slope 
of —23.3 ± 3 s~^. (b) Average on-time of events whose trans- 
lations fall within the same narrow range. The inset presents 
the same data on a In-linear scale, the solid line being an 
error-weighted linear fit with a slope of 0.040 ± 0.004 s. 



value of n, Pne{n,e) = [n/(27rcr^i)]i/2e-"eV(2<T,\), This 
implies that, regardless of the distribution of the aver- 
age absolute angle of events having the same translation 
I should decrease as 

\ei\ = [2al,h/{nl)Y/\ (1) 

which is consistent vi^ith the findings of Fig. Ol^b) . From 
the fit we get a^^h = 500 ± 200 (°)2/im. Thus, large 
cascades give rise to "focusing" onto small-angle trans- 
lations. This effect, along with the broad I distribution, 
explains the large statistical weight of small angles. 

Suppose that each fold takes a time Ti to complete 
and another time r to "topple" another fold. Consider 
a cascade made of g generations of topples, where each 
fold can topple q others |0. The total on-time is T = 
Ti + {g — 1)t, and the total number of folds in n — g ii 
q = 1, or = (qf — l)/(g — 1) if g > 1. From these two 
equations and n ^ l/li we obtain a relation between the 
on-time T and the translation I, T = {t/Ii)1 + Ti — t 
if (? = 1, and T = (r/ Ing) ln[((j - + 1] + Ti - t if 

q > 1. Thus, for low cooperativity [q — 1) we expect T to 
depend linearly on I, whereas for larger cooperativity the 
dependence becomes logarithmic. As seen in Fig. 0Jb), 



the results are consistent with the latter, i.e., the cascades 
are cooperative, one fold toppling several others. 

Assuming a fixed probability a to topple a genera- 
tion of folds, we get an exponential distribution of g, 
Pg{g) = [(1 — a)/a\a^ , and hence an exponential distri- 
bution of T for T > Ti. The statistics for T < Ti are 
determined by the scatter of Ti , which has been neglected 
so far. Hence, pt{T) should have an asymmetric shape, 
dropping sharply for T < Ti and decaying exponentially 
for T > Ti. These conclusions agree with the measure- 
ments of Fig. 01a), whereupon the distribution peak can 
be identified as the single-fold time, Ti = 0.12 ± 0.03 
s. From the first two moments of pt{T > Ti) we ex- 
tract T = 0.026 ± 0.005 s and a = 0.54 ± 0.05. Hence, 
{g) — 1/(1 — a) = 2.2 ± 0.2, i.e., the cascades consist of 
2-3 generations on average. From the linear fit in Fig. 
^b), whose slope is equal to t/ In q, wc get q ~ 1.9 ± 0.2, 
i.e., each fold topples about two others. The cooperativ- 
ity and short toppling time account for the uniformity of 
event duration — even the largest cascades involve only a 
few generations and do not last long. 

Using the exponential distribution of g and the relation 
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between g and n, we find 

p„(n) ~ [(q- l)n+l]-'5, /? = 1 - In a/ In g. (2) 

The distribution of magnitudes tends for large n — l/li 
to a power law, the decay exponent being always smaller 
than —1. (For our system we get (3 = 2.0 ± 0.2.) Using 
the expressions for pn{n) and Pn0{Q) we can calculate the 
distribution of angles as = 'Y^°^=iPn{''^)Pne{Q), yielding 

oo 

pom) = ^0 - 1)" + l]-^n''^e--'^'^^<\ (3) 

n=l 

where Me = 2{2nal^)-^/^{q ~ lf/aP,q/{q - 1)] is a 
normalization factor. This distribution reproduces well 
the focusing of the angular distribution, as shown by the 
solid line in Fig. Eta). We note that, if one has in a 
certain system 1 < < 3/2 (i.e.. In a/ In q > —1/2), Eq. 
© predicts an (integrable) singularity of at 6* — > 0. 
The fit in Fig. |2Ia) gives aei — 16 ± 2° which, together 
with the estimate for cr'^ili obtained from Fig.l^Ib), yields 
li — 2.0 ± 0.8 /^m. Finally, we use the values derived for 
li, q and (3 to reproduce the translation distribution pi 
according to Eq. Q. As shown by the solid line in Fig. 
|2Ia), the calculated distribution gives a reasonable fit to 
the decaying part of pi{l). (The sharply increasing part 
is probably determined by the scatter of q.) 

The cascade analysis provides a consistent account of 
the measured statistics. The alternative single-fold sce- 
nario, corresponding to n = g = 1, does not agree with 
the measurements and cannot account for the sharp dis- 
tributions of angle and duration. The evidence for co- 
operative cascades, however, remains indirect. Our anal- 
ysis implies that the folding transition follows unusual 
nucleation kinetics, in which single-fold growth is macro- 
scopic in one dimension (length) but restricted in another 
(width) . Consequently, a single nucleus cannot fully relax 
its super-stressed environment, thereby driving the nu- 
cleation of other folds in a chain-reaction manner. Such 
a process, resembling an autocatalytic chemical reaction 
or nuclear fission, has never been recognized, to the best 
of our knowledge, in the context of mechanical instabil- 
ity such as the buckling discussed here. Moreover, this 
scenario should not be restricted to our specific system 
but is to be expected whenever there is an autocatalytic 
instability whose evolution is limited for some reason to 
discrete units of relaxation. 

A key question is what sets the scale of the restricted 
fold growth. Folding may introduce extra strain in the 
monolayer, e.g., as a result of mismatches in the rapidly 
folded region. This strain will increase with fold width, 
eventually balancing the cohesion energy and halting the 
folding. Another open issue is the mechanism of correla- 
tion between folds in a cascade. Toppling may be caused 
by the long-range stress field emanating from the tips of 
a propagating fold. This extra stress appears immedi- 
ately after fold nucleation and can account for the short 



toppling time inferred above. A complicated question to 
be addressed in a future publication j^j relates to the role 
of compression rate. It will be interesting to check what 
happens to the folding cascades when the compression is 
not unidirectional, e.g., in a circular trough |l2| . 
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